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We present results for the renormalised light and strange quark masses calculated using Domain Wall Fermions 
in quenched QCD. New results using the DBW2 gauge action at « 2GeV and ~ l.SGeV will be presented 
and compared against existing results at ~ 2GeV using the Wilson gauge action. This comparison allows a 
study of the uncertainties due to both finite lattice spacing and residual chiral symmetry breaking effects. 



1. INTRODUCTION 

We report here results for the light and strange 
quark masses extracted from the quark mass de- 
pendence of the meson spectrum calculated us- 
ing domain wall fermions [Qj^ (DWF) , see ^ for 
our notation and conventions, in the quenched 
approximation and renormalised using the non- 
perturbative renormalisation (NPR) technique 
developed by the Rome-Southampton group 

All the results presented will be from sim- 
ulation using a lattice with spatial dimensions 
16'^ X 32 and a fifth dimension of size 16. The new 
results use the DBW2 gauge action [|| at /3 = 0.87 
and (3 = 1.04 corresponding to a"^ « 1.3GeV 
and a^^ « 2GeV respectively. This action is cho- 
sen for its much smaller explicit chiral symme- 
try breaking when compared to other actions 
The results will be compared with those of pre- 
vious studies using the Wilson gauge action with 
P = 6.0 §,0j|. For P = 0.87, M5 = 1.8 was 
used and 53 and 100 configurations were collected 
for the NPR and hadron spectrum respectively, 
whereas for /3 — 1.04, M5 = 1.7 was used and 51 
and 405 configurations were collected. 

2. RENORMALISATION 

Our final result will be the renormalised mass 
in the MS'-scheme at 2 GeV. The relation be- 
tween this value and the bare mass, uif, consists 



of several factors and may be written 



= Z 



match r^BJ 
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The residual mass, TOrcs, occurs because of small 
explicit chiral symmetry breaking effects due to 
the finite extent of the fifth dimension. This we 
determined from the breaking term of the axial 
Ward-Takahashi identity |^. is the mass 

renormalisation factor in the RI/MOM-scheme 
which is calculated directly on the lattice using 
the method of the Rome-Southampton group jj] 
and then matched to the M5-scheme at 2 GeV 
using a two-loop perturbative matching factor, 

2-match 

For the NPR technique to be successful a "win- 
dow" of momenta must exist for which contam- 
ination due to (low-momenta) condensate effects 
are small but also for which (high-momenta) 
lattice artifacts are suppressed. Previous work 
shows that when using DWF with the Wilson ac- 
tion at a^^ « 2GeV this window exists |8j. In this 
work we confirmed that this is also the case for 
the DBW2 action and checked the feasibility of 
NPR at a-^ « 1.3GeV. 

While Zm may be calculated directly from 
the momentum-space quark propagator this ap- 
proach has large systematic errors. Instead we 
consider the renormalisation of the flavour non- 
singlet local quark bilinear operator, uTd. The 
NPR technique provides a clean extraction of the 
ratio Zr/Zq, where Zq is the quark renormalisa- 
tion factor. Putting the value of Zs/Zq (scalar) 
and ZA/Zq (axial) together with the value of Za 
calculated from the hadronic matrix element of 
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Figure 1. ZsjZ^ bare and RGI for the DBW2 
action at /? = 1.04 



Figure 2. ZsjZ^ bare and RGI for the DBW2 
action sX (i — 0.87 



the focal and conserved axial currents with the 
pseudo-scalar density]^ gives our final value for 
Zg. This can then be used to calculate Z^ as 
Zs = l/Z,^. 

An important check of the existence of a "win- 
dow" may be made by dividing ZgjZq by the 
predicted running from 3-loop perturbation the- 
ory to produce a renormalisation group invari- 
ant (RGI) value. Fig. || shows this for /3 = 1.04 
DBW2 where the RGI and bare data are con- 
strained to agree at (ap)^ — 1. 

Fig. I show this for ^ = 0.87 DBW2 (in this 
case the bare and RGI are constrained to agree at 
(ap)^ = 1.5). There is both more evidence of cur- 
vature at low-momenta due to condensate effects 
and a noticeable slope versus (ap)^ at high mo- 
menta. Interpreting this slope as O(a^) contam- 
ination and fitting the RGI value to B{apf 
allows the final value for Zs/Zq to be extracted. 

Table |^ collects together the (preliminary) 
DBW2 results for Za and Zs/Za (the latter in 
the MS-scheme at 2 GeV), and the results for 
Wilson gauge results already quoted in j^. 



Table 1 



Results for Za and Zs/Za 



Action 


P 


Za 


Zs/Za 


DBW2 


0.87 


0.7776(5) 


0.898(16)(90) 


DBW2 


1.04 


0.8402(2) 


0.849(15)(40) 


Wilson 


6.00 


0.7555(3) 


0.830(09)(35) 



3. MESON MASSES 

For the pseudo-scalar meson masses chiral per- 
turbation theory to first order tells us that 

m^a^ = -B^am (2) 
TO^a^ = BTra{ms +rn)/2 (3) 

where m — ^{mu+md), while for the vector 
mesons the leading dependence should be linear 
in the quark masses 

iripa — Ap + BpYn . (4) 

Our approach to calculating the quark masses is 
to first fit the meson spectrum for the coefficients 
Stt, Ap and Bp and then to use the physical values 
for m,r, and rup to calculate a, rn and mg. 

In extracting the meson masses we must deal 
with two systematic effects due to working in 
the quenched approximation: the contamination 
of meson correlators due to the presence of un- 
suppressed topological near zero- modes and 
significant deviations from Eq. ^ at light masses 
due to quenching effects. 

The effects of zero-modes are most apparent 
at light masses and small volumes, but also are 
strongly dependent on which correlators are used 
to calculate the mass. In particular when ex- 
tracting the pseudo-scalar meson mass the zero- 
mode effects may be reduced (but not completely 
avoided) by calculating the mass from axial cor- 
relators rather than the more usual pseudo-scalar 
correlators. In addition to using the axial cor- 
relator for the pseudo-scalar mass extraction, we 
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Figure 3. Pseudo-scalar mass versus m/ + 
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Figure 4. Pseudo-scalar mass versus renormalised 
mass 

only use bare masses greater than nif = 0.01 in 
our analysis. 

Using the mass extracted from the axial cor- 
relator, the consistency of our different data-sets 
with each other and the NPR calculation may 
be tested. Fig. ^ shows the pseudo-scalar mass 
versus rh — mj + rrircs, both measured in MeV, 
whereas Fig. ^ shows the same data versus renor- 
malised mass. As can be seen once the renormali- 
sation factors are taken into account the data-sets 
show excellent consistency both between the two 
gauge action at the same lattice spacing and the 
two different lattice spacing with the same gauge 
action. The deviation from Eq. ^ is very appar- 
ent in our data with a simple linear fit of 



(5) 



for the (3 — 1.04 DBW2 data having a constant 
term, inconsistent with zero by 3 standard 



deviation whereas the intercept of a more general 
quadratic fit misses zero by 7 standard deviations. 

As this is the case we include a quenched chiral 
logarithm in our fitting formula (but drop its ef- 
fects for the final calculation of the quark masses). 



B-jrfn 1 — (51oe 



~A2- 



(6) 



where we use a chiral scale of A = IGeV [|lO| . 
This gives very consistent values of (5 = 5.1(19) x 
10^2 and 5 = 4.6(18) x 10"^ foj. Wilson and 
DBW2 at 2GeV and 5 = 2.28(12) x 10"^ at 
1.3GeV. However, none of the fits are of high 
quality and the best treatment of the quenching 
effects is an issue that requires further study. 

4. RESULTS 

Putting together the NPR results and the fits 
to the meson spectrum gives values for the quark 
masses in MS'-scheme at 2 GeV. Using Equations 
to ^ implies that 

m,/m 2m\/ml - 1 « 26 (7) 
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exactly for our results so we only quote the 
strange quark mass. For this we get 132.5(31) 
MeV (/? = 0.87, DBW2), 132.6(21) MeV {(3 = 
1.04, DBW2), and m, = 125.6(28) MeV (/3 = 6.0, 
Wilson) where the errors quoted are purely sta- 
tistical. 
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